Exact master equations for the non-Markovian decay of a qubit 
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Exact master equations describing the decay of a two-state system into a structured reservoir are 
constructed. Employing the exact solution for the model we determine analytical expressions for the 
memory kernel of the Nakajima-Zwanzig master equation and for the generator of the corresponding 
time-convolutionless master equation. This approach allows a detailed investigation and compari- 
son of the convergence behavior of the corresponding perturbation expansions. Moreover, we find 
that the structure of widely used phenomenological master equations with memory kernel may be 
incompatible with a non-perturbative treatment of the underlying microscopic model. We discuss 
several physical implications of our results on the microscopic analysis and the phenomenological 
modelling of non-Markovian quantum dynamics of open systems. 
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I. INTRODUCTION 

The field of open quantum systems is of great 

interest because of its relevance in applications of quan- 
tum mechanics, as well as in a deeper understanding of 
the theory itself. Indeed the study of the interaction 
between a quantum system and its environment is an 
endeavor common to many fields such as quantum mea- 
surement theory, quantum communication, quantum op- 
tics, condensed matter theory and quantum chemistry to 
name a few. The field is well assessed as far as Markovian 
dynamics is concerned, in which the Gorini-Kossakowski- 
Sudarshan-Lindblad expression for the generator of a 
quantum dynamical semigroup 0, H| provides a bench- 
mark result for both microscopic and phenomenological 
approaches. This situation is however not satisfactory 
when one has to go beyond the Born-Markov approxima- 
tion and considers systems in which a separation of time 
scales between system and environment can no longer be 
assumed in a realistic description. Memory effects then 
become important and a non-Markovian description is 
mandatory. For this case a general consistent theoreti- 
cal framework has not yet been found, and partial results 
have been obtained as a result of intense efforts (see, e.g., 
Refs. An important step in the development of a 

general theory consists in the construction of a suitable 
measure that quantifies the degree of non-Markovianity 
for a given dynamical evolution [22, [23| . 

In this article we will obtain the exact Nakajima- 
Zwanzig kernel for a two-level system coupled to a 
Bosonic reservoir discussed in 0], and compare it to the 
exact time-convolutionless master equation as well as to 
the Markovian approximation of the dynamics. This 
will show how involved the transition from the approx- 
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imate Markovian level of description to the exact non- 
Markovian regime can be. Indeed, the non-Markovian 
memory kernel is found to have an operator structure 
which differs from the one that appears in the Born- 
Markov approximation. Often one tries to obtain dynam- 
ical equations of motion for non-Markovian systems by 
slight modifications with respect to the Markovian case, 
e.g., by considering a master equation which involves a 
superoperator given by a convolution in time of the corre- 
sponding Markovian superoperator (2^ - l29| . Our results 
show that such an approach, although being justified as 
a phenomenological modelling, can be incompatible with 
a non-perturbative treatment of the underlying micro- 
scopic system-environment model. Moreover, different 
perturbation expansions such as time-convolutionless and 
Nakajima-Zwanzig projection operator technique turn 
out to have different ranges of validity. Indeed the time- 
convolutionless expansion breaks down at finite time in 
the strong coupling limit, while the Nakajima-Zwanzig 
approach does not preserve positivity if restricted to sec- 
ond order. Furthermore the convergence to the exact 
solution is not uniform with respect to the expansion pa- 
rameter: Different matrix elements of the statistical op- 
erator such as coherences and populations are obtained 
with quite different accuracy at the same perturbative 
order. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the model and its exact solution, which is later 
exploited to obtain the exact equations of motion for the 
reduced statistical operator of the system. In Sec. IIIII 
we recall the structure of the time-convolutionless mas- 
ter equation, pointing out two different perturbation ex- 
pansions for the generator. In Sec. IIVI we derive the 
Nakajima-Zwanzig integral kernel, providing an alterna- 
tive expansion with respect to the standard method. The 
two results are compared in Sec. El also building on an 
exact analytic expression for all the quantities involved 
obtained considering a Lorentzian spectral density. We 
finally draw our conclusions in Sec. I VI I 
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II. THE MODEL AND ITS EXACT SOLUTION 

The total Hamiltonian of the model is given by 

H = Hs + H E + Hi = H + H I , (1) 

where 



H s = loq<j + (7- 



(2) 



describes a two-state system (qubit) with ground state 
0), excited state |1) and transition frequency ujq. The 
operators er + = |1}(0| and <r_ = |0)(1| are the raising 
and lowering operators of the qubit. The environmental 
Hamiltonian is taken to be 



H E = ^uj k blbk, 



(3) 



describing a collection of harmonic oscillators with cre- 
ation and annihilation operators b\ and bk which satisfy 

bosonic commutation relations [bk,b\.,] = Skk' ■ The in- 
teraction Hamiltonian takes the form 



Hj = ^ (gk(T+ <8> b k + g* k o_ <g> b\ 



(4) 



The model thus describes for example the coupling of 
the qubit to a reservoir of electromagnetic field modes 
labelled by the index k with corresponding frequencies 
uj k and coupling constants g k , and has already been dis- 
cussed in 

In the following we will work in the interaction picture 
with respect to Hq = Hs + He- In this picture the 
Schrodinger equation reads 



-|*(t)> = -ifl,(f)|*(t)>, 



(5) 



where the interaction Hamiltonian is given by 

H I {t)=a + {t)®B(t)+<7-(t)®B' ! (t) (6) 

with 

a ± (t) = a ± e ± ^ t (7) 

and 

B{t) = Y,9kb k e- i " ht - (8) 



It is easy to verify that the operator 



N = a+a- +XX 6fc 



(9) 



for the number of excitations in the system commutes 
both with the total Hamiltonian H and with the inter- 
action Hamiltonian Hi(t), which is a consequence of the 
fact that the rotating wave approximation has been used 



in the interaction Hamiltonian ((4]). It follows that any 
initial state of the form 



|tf(0)) = co|0>® |0}jj + ci(0)|1)®|0}js 
+ 5> fc (0)|0>® |fc) E 



(10) 



evolves after time t into the state 



|*(t)> = Co|0}® |0>B + Cl(f)|l)®|0} B 

+ J2 c k(t)\0)®\k) E - 



(11) 



The state \0}e denotes the vacuum state of the reservoir, 
and \H)e — &1|0}_e the state with one particle in mode 
k. Note that the amplitudes C\{t) and c k (t) depend on 
time, while the amplitude cq is constant in time because 
of i?j(*)|0) ® |0)_b = 0. Substituting Eq. JlTJ into the 
Schrodinger equation <j5j) one finds 



T Mt) = 



dt 
d_ 



9ke 



i(wo— Wfe)< 



Ck(t), 



= —i n* p-i^o—^kft 



c k (t) = -ig* k 



Cl(t). 



(12) 
(13) 



We assume in the following that c/-(0) = 0. This means 
that the environment is in the vacuum state initially and 
that the total initial state is given by the product state 

|*(0)) = (c |0) + ci(0)|l)) ® \0) E = IV'(O)) ® \0)e- (14) 



Expressing c k {t) in terms of C\(i) by means of Eq. {Te 
and substituting the result into Eq. (fT2"j) one obtains an 
integrodifferential equation for the amplitude Ci(i), 



d , , 
dt C ^ 



dtif{t-tx)ci{ti). 



(15) 



Given the solution of this equation, which can be found 
through a Laplace transformation, the amplitudes Cfc(i) 
are determined by Eq. (IT31 . The kernel fit — t\) of 
Eq. (|15|) is given by a certain two-point correlation func- 
tion of the reservoir, 

f(t-h) = {0\B(t)B^t 1 )\0)e^ o(t - tl '> (16) 



k 



| 5fe |2 e i(u )0 - w , ! )(t-t 1 ) ) 



on which no restrictive hypothesis is made, so that our 
results will be valid for an environment with a generic 
spectral density. 

With the help of the procedure described above, al- 
ready used by Weisskopf and Wigner in their classical 
paper (531 , one finds the solution of the Schrodinger equa- 
tion of the total system with initial states of the form 
(TlH) lying in the sector of the Hilbert space correspond- 
ing to zero or one excitations. By means of this solution 
we can construct the exact dynamical map describing the 
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time-evolution of the reduced density matrix of the qubit 
which is given by 

*)=»«w»=(S;li)- ("> 

where Pij(t) = (i\p(t)\j) for i, j = 0, 1. Using Eq. (|TT|) we 
find 



pu(t) = l-poo(t) = |ci(t)| 2 , 
Pio(*) = = cSci(i). 



(18) 
(19) 



It is convenient to define the function G(t) as the solution 
of the equation 



tfti/(t-ti)G(ti) 



(20) 



corresponding to the initial condition G(0) = 1. We then 
have c\ (t) — G(t)c\ (0) and, hence, the dynamics of the el- 
ements of the reduced density matrix can be represented 
as follows, 

p u (i) = \G(t)\ 2 Pll (0), (21) 
Poo(i) = poo(0) + (l-|G(i)| 2 )pii(0), (22) 
fHo(t) = G(f)pio(0), (23) 
PoiW = G*(t)p 01 (0). (24) 

These equations have been derived for the pure product 
initial state (fl"4")l . i.e., they describe the time-evolution 
corresponding to the pure reduced system's initial state 
p(0) — \ip(0))(4>(0)\. However, since any mixed initial 
state can be represented as convex-linear combination 
of pure initial states, and since the function G(t) intro- 
duced above docs not depend on the initial condition, the 
equations (|2ip - (|24[) hold true for any pure or mixed ini- 
tial state. They thus represent the exact dynamical map 
&(t) which transforms the initial states into the states at 
time t: 



*(t) : p(0) p(t) = *(f)p(0), 



t > 0. 



(25) 



Since we have constructed this map from the exact solu- 
tion of the model, it is clear from the general theory of 
open quantum systems that is completely positive 
and trace preserving. 



III. THE TIME-CONVOLUTIONLESS MASTER 
EQUATION 

A. Exact master equation in time-convolutionless 
form 



governing the dynamics of the reduced density matrix. 
The time-convolutionless generator is defined in terms of 
the dynamical map $(i) by means of 



/C TCL (t) = Ht)$- l {t) 



(27) 



provided the inverse map $ _1 (i) exists. Using then 
Eqs. (f2~T j) - (|2"4")) one shows that the generator takes the 
following form [2j, 



£tCl(*)p = —=S{t)[<J + a_,p} 



(28) 



-lit) 



a -pa ^ 



{a + a-,p} 



where we have introduced the definitions 



lit) 



-25ft 



G(f) 
Git) 



Sit) 



-23 



Ojt) 
Git) 



(29) 



By construction, Eq. (|2"6]) with the generator (|2"8")l repre- 
sents an exact time-local master equation. Note that the 
generator is well-defined as long as G(t) ^ 0. The quan- 
tity S(t) plays the role of a time-dependent frequency 
shift, and j(t) can be interpreted as a time-dependent 
decay rate. We observe that the structure of /Ctcl is 
similar to that of a Lindblad generator. However, due 
to the time dependence of the coefficients S(t) and j(t) 
Eq. (|26[) does generally not yield a quantum dynami- 
cal semigroup. Moreover, the time-dependent rate "f(t) 
may become negative, signifying strong non-Markovian 
behavior of the reduced system dynamics. 



B. Perturbation expansions of the generator 

In most cases of interest the time-convolutionless gen- 
erator can only be determined through a perturbation 
expansion. Here we investigate two methods of expand- 
ing the exact master equation (|2"oT) with respect to the 
strength of the interaction Hamiltonian Hj. To this 
end, we introduce a small overall expansion parameter 
a, replacing the coupling constants in the interaction 
Hamiltonian Q by agk- The two-point correlation func- 
tion f(t), being proportional to a 2 , is then to be regarded 
as a quantity of second order. 

The first method consists in using Eq. (I20I) to obtain a 
perturbative expression for G(t) from which one directly 
finds an expansion for the coefficients j(t) and S(t) ap- 
pearing in the master equation. The expansion of G(t) 
is obviously of the form 



G(t) = J2a 2n G (2n) (t), 



(30) 



The exact solution determined in Sec. |TT] enables the 
construction of the exact generator /Ctcl of the time- 
convolutionless master equation 



dt 



p{t) = /Ctcl (*)p(t) 



(26) 



where G^°'(t) = 1 because of the required initial condi- 
tion G(0) = 1, and Eq. (|2H)) leads to the following recur- 
sion relation 



G (2n) (t) = - 



t rtl 

dh / dt 2 f(h 
ii Jo 



i 2 )G (2 "~ 2) (i 2 ). (31) 
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To illustrate the procedure we determine the frequency 
shift and the decay rate to fourth order in a: 



(32) 



= a 2 G^(t)+a 4 [G^(t)-G i2 Ht)G i2 Ht)\ +0(a 6 ). 

With the help of these expressions one obtains the second 
and the fourth order contributions for the coefficients of 
the master equation: 



(33) 



7 (4) (i)+iS (4) (<) = 2 f dtt f 1 dt 2 f 2 dt 3 
Jo Jo Jo 

x [f{t - t 2 )f(h - t 3 ) + f(t - t 3 )/(ti - t 2 )\ . 

Another possibility for the construction of the pertur- 
bation expansion is to use the general method of ex- 
panding the time-convolutionless generator in terms of 
the ordered cumulants. This procedure allows to write a 
closed expression for the coefficients of the master equa- 
tion which takes the form (for details, see and refer- 
ences therein): 

7 (2n) (t) + iS ( - 2n \t) = [ dh ^ dt 2 ...[ 2n 2 dt^-i (34) 
Jo Jo Jo 

X2(~l) n+1 (f(t - h)f(t 2 - t 3 ) . . . /(t 2 „_ 2 - * 2 „_l)) oc . 



IV. THE NAKAJIMA-ZWANZIG MASTER 
EQUATION 

A. The exact memory kernel 

The Nakajima-Zwanzig master equation is given by 



dt' 



p(t) = / dtiK. N z{t -ti)p(ti), 



(35) 



where the superoperator /Cnz(t) represents the memory 
kernel. We construct the form of this kernel from the 
exact solution of our model obtained in Sec. HU To this 
end, we employ the following Ansatz, 



/Cnz(t)p = -ie(T)[(j + (T-,p\ 
+k 1 (r 



(36) 



-k 2 (r) 



1 r 

cr+a-pa+a- - - {(T + (T- 



where the functions e(r), fci(r) and k 2 (r) are real, such 
that the master equation preserves Hermiticity and trace. 
The equations of motion for the population pu(t) and 
the coherence pw(t) obtained from this master equation 
read: 



(37) 



and 

d 
dt 



Pw(t) = 



dti 



-{k^t - h) + k 2 (t -h)} + ie(t - ti) 



(38) 



On the other hand, Eq. Q together with Eq. (J20]) yields: 



dt 



Pio(t) 



dt 1 f(t - ti)p w (ti), 



(39) 



where we have set the expansion parameter a equal to 
one and we only have to remember that f(t) is a quantity 
of second order. Comparing Eqn. (I3"i?l) with Eq. we 
see that the expression within the square brackets of (|38|) 
must be equal to f(t — t\), i. e., we get the conditions: 



e{r) 

fci(r) + fc 2 (r) 



/ 2 (r), 
2/i(r), 



(40) 
(41) 



where /i(r) and f 2 (r) denote the real and the imaginary 
part of the correlation function: 



/(r) = /i(r)+i/ 2 (r). 



(42) 



In order for Eq^. (1371) to reproduce the correct solution 
([2~T|) we have to choose k\{r) in such a way that the so- 
lution of the equation 



d 



z(t) 



dt ' 
is given by 



dt\ki{t — ii)z(ti), 



z(0) 



z(t) = \G{t)\ 2 . 
Formulated in Laplace space this means that 

I — uz(u) 



ki (u) 



z(u) 



(43) 



(44) 



(45) 



Since the superoperator (|36[) preserves the Hermiticity 
and the trace of the density matrix, Eqs. (|2II) - (f2~4")) fol- 
low directly from Eqs. (|37|) and (|38|) . Thus, we find 
that Ec|. p6p represents the exact memory kernel of the 
model for any given two-point correlation function. In 
fact, given /(r), the functions e(r), fci(r) and k 2 (r) are 
uniquely determined by Eqs. (|40)) , (I4T1) and (|45")) . In view 
of this result the memory kernel (|36|) can now be written 
in the form 

1Cnz(t)p = -if 2 {r)[a + a_,p\- fi{r){a + a_,p} 

+ ki(T)(J_p<7 + 

+ [2f 1 (r)-k 1 (T)]a + a_pa + a_, (46) 

which only involves the real and the imaginary part of the 
correlation function and the function ki (r) which has to 
be determined from Eq. (|45[) . 
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We note that the various coefficients in the memory 
kernel exhibit a very different convergence behavior. In 
fact, we see that the commutator and the anticommuta- 
tor term in Eq. (|46[) come out exactly in second order 
in a. It follows that the equation of motion for the co- 
herence pio [see Eq. ([39]) ] is already reproduced exactly 
within second order, while the exact representation of the 
equation for the population pn requires in general the in- 
clusion of all orders of the expansion. This non-uniform 
convergence behavior of the elements of the density ma- 
trix has been observed also in other, more complicated 
models [3l[ , and seems to be a typical feature of the per- 
turbation expansion of the memory kernel. 

As will be shown below the relations (|20|) together with 
(|43| - (j45| provide a direct perturbation approach to the 
determination of the functions appearing in the mem- 
ory kernel Eq. (l36l) . as an alternative to the standard 
Nakajima-Zwanzig perturbation expansion. Moreover, 
this set of equations allows in some cases to derive a 
closed analytical expression for the memory kernel. 



B. Perturbation expansions of the memory kernel 

Here we discuss two methods of expanding the exact 
memory kernel with respect to the strength of the interac- 
tion Hamiltonian Hi. The first expansion method relies 
on the expansion Eq. (1501) for the solution of Eq. (1231) . 
which determines the dynamical map once the two- 
point correlation function f(t) of the model given by 
Eq. (flB) is specified. 

Indeed as shown in Sec II V Al to obtain the memory ker- 
nel Eq. (|33|) we only need to know the function ki(t). A 
perturbative expression for the latter can be easily ob- 
tained relying on the expansion Eq. (|30|) for G(t), noting 
that thanks to Eq. (1431) the Laplace transform of k% (t) can 
be directly expressed through Eq. (1431) by means of the 
Laplace transform of the function z(t) = |G(t)| 2 . This 
procedure leads to the following expansion 



hit) 



n=0 



(47) 



as described in detail in Appendix [Al where the zero 
order contribution is immediately seen to be zero. 

Here we consider for the sake of simplicity only the 
second order contribution, which is readily obtained. Ac- 
cording to Eq. (|3"Tj) together with the initial condition 
G(0) = 1 the expression for G(t) up to second order is 
given by 



G(t) 



1 



t rti 

dh / dt 2 f{t 2 ), 
o Jo 



(48) 



so that in the same approximation, recalling that the 
two-point correlation function fit) is a quantity of second 
order, one has 



(49) 



*(t)«l-2 / dti / dt 2 /i(*a) 



where according to Eq. (|42j) f\ it) denotes the real part of 
the correlation function. The Laplace transform of this 
quantity is now easily expressed in terms of the Laplace 
transform of the correlation function according to 



z{u) 



«-2/i(u) 



and further exploiting Eq. (|45[) we find 
fci(tt)«2/i(«). 



(50) 



(51) 



This immediately implies for the second-order contribu- 
tions to the kernel (1331) 



(52) 



^ 2) (i) = 2/ 1 (r), 
and therefore due to Eq. (j4T]) 

fcf (r) = 0. 



(53) 



As shown in Appendix [A] the fourth-order contribution 
reads 

kf\t-h) = -25R f dt 2 f 2 dt 3 [fit - t 3 )f(h - h) 
+ f{t-h)fih-t 2 )], (54) 



so that = 



k[ 4 \ Indeed Eq. (|41[) generally implies 

k^ n ^ — —k^' for n > 2, therefore Eq. (l47l) also provides 
an expansion for k 2 it). 

The second expansion method is to employ the gen eral 
Nakajima-Zwanzig projection operator technique [32I.I33T] 
in which the memory kernel is expressed in terms of 
the full propagator of the total system. The details of 
this method for our model are presented in Appendix [Bj 
where it is shown that the projection operator technique 
reproduces, as expected, the above results obtained by 
the direct expansion of the coefficients in the memory 
kernel. 



V. DISCUSSION 

A. Comparison of the time-convolutionless and the 
Nakajima-Zwanzig master equation 

It is interesting to compare the time-convolutionless 
master equation Eq. (1261) with the Nakajima-Zwanzig 
master equation Eq. (|35p . For the considered model the 
functions appearing in Eq. (|28p and Eq. (|36p are given 
by Eq. (|2S} and Eqs. ([4Tj |) -([4T )l respectively. 

Comparing Eqs. (|28p and (133)) we see that the super- 
operator structure of the memory kernel differs from that 
of the time-convolutionless generator. In fact, the mem- 
ory kernel (|36p contains the term proportional to k 2 ir) 
which involves the projection (7 + cr_ = |1)(1| onto the 
excited state. Without such a term the equations (|3"?| 
and ([33j) for the population and the coherence would be 
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incompatible with the exact expressions (|2Tj) and (1231) . 
However, a term with this structure is missing in the 
time-convolutionless generator ([28)) . A further remark- 
able point is the fact that in second order /^(t) = 0, 
according to Eqs. (|4Tj) and (f52j) . This shows that the 
difference in the superoperator structure of the memory 
kernel and the time-convolutionless generator is visible 
only in higher orders of the perturbation expansion. 

The above discussion leads to some conclusions which 
are important for the modelling of non-Markovian dy- 
namics through phenomenological master equations. In 
the Markovian limit our model yields the following Lind- 
blad generator C describing a quantum dynamical semi- 
group, 



Cp 



-Sm{<J+<J-,p\ + 7m 



a- pa- 



,P} 



(55) 

with constant frequency shift Sm and decay rate jm > 0. 
Usually master equations of this form are derived by ap- 
plying the Markov approximation and second order per- 
turbation theory (Born-Markov approximation). A natu- 
ral and widely-used non-Markovian generalization is then 
obtained from this equation by keeping the structure of 
the Lindblad generator £ and by introducing a certain 
kernel function h(r) to arrive at a master equation of the 
form iJi 



d , . 

5' w 



dhh(t-ti)Cp(ti). 



(56) 



Although this equation is perfectly justified as a phe- 
nomenological ansatz, it does in general not represent 
the correct structure of the memory kernel of the under- 
lying microscopic model. In fact, we see that even for the 
simple model studied here the true memory kernel (|36j) is 
not of the form of Eq. (|56[) . but involves additional terms 
that are absent in the Markovian Lindblad generator. In- 
deed it is rather given by a linear combination of terms of 
the form Eq. (|56[) , where going beyond the Born approx- 
imation besides the Markovian Lindblad generator other 
operator structures appear, which are still in Lindblad 
form but with different Lindblad operators. This obser- 
vation seems to be of particular relevance for the analysis 
of the positivity and the complete positivity of the dy- 
namical maps obtained from phenomenological equations 
of motion. 



B. Example 

These considerations can be nicely illustrated consid- 
ering the example of an exponential correlation function, 
corresponding to a Lorentzian spectral density [2j 



fir) 



(57) 



where the parameters 70 and A are real and positive. 
For this case both time-convolutionless generator and 



Nakajima-Zwanzig kernel can be exactly calculated. In- 
deed by means of Eq. (|20p one obtains for the function 
Git) the expression: 



G(t) = e 



- P -M/2 



cosh ( -^-fij + h sinh ( ^-8 



(58) 



where 8 = yl — 2jq/X. Note that this function is always 
real. Furthermore it stays positive for any time t in the 
weak coupling regime 70 < A/ 2, while for strong coupling 
70 > A/2 the parameter 8 becomes purely imaginary and 
the function G(i) starts to oscillate. In particular it goes 
through zero for the first time when t is equal to the 
smallest positive solution of 



2 

aI 



arctan 8 



(59) 



where 8 = ^270/ A - 1 and ntN. Building on Eq. (|58|) 
one can obtain the exact expressions for the functions 
j(t) and S(t) appearing in the time-convolutionless gen- 
erator, given by S(t) = and 



sinh (%6) 



,5 cosh (f 8) + sinh(f 8) 



(60) 



In order to obtain the Nakajima-Zwanzig kernel one 
considers the Laplace transform of the function z(t) = 
|G(t)| 2 which is found to be: 



z(u) 



(u+ A)(u + 2A) +70 A 



(61) 



(u + A) [(« + A) 2 - A 2 + 2 7o A] ' 
so that according to Eq. (1151) one has: 

fci(«) = 7bA 7 ^ "^u x - ( 62 ) 
[u + A)(u + 2A) + 70 A 

Transforming back to the time domain we finally get 



fc lW = 70 A e - 3 ^ 2 



cosh 



Af A 1 . , f\t., 
8 I + — sinh I —6' 



2 ~ J ' 8' ~~~ V 2 

' (63) 

where 8' = yl — 470/A. Substituting this result into 
Eq. ()46|) we find the exact memory kernel for the case of 
an exponential correlation function. 

The exact expressions Eq. and Eq. (153")) already 
allow for an important comparison. While the func- 
tion on the right-hand side of Eq. (j63p represents an 
analytic function of 70 (remember that 70 is a quan- 
tity of second order in the expansion parameter a), so 
that the Nakajima-Zwanzig memory kernel has an infi- 
nite radius of convergence, the same does not hold true 
for the time-convolutionless generator. Indeed the time- 
convolutionless expansion breaks down in the strong cou- 
pling regime 70 > A/2 when the function G{t) given in 
Eq. (j58[) goes through zero, corresponding to the diver- 
gence of the decay rate 7(f) given in Eq. (jSO 7 )) and ob- 
tained from the relation Eq. (l2l)|) . 
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Considering an expansion in 70 of the function ki(t) 
given by Eq. (I63[) which fixes the memory kernel, due to 
the fact that the correlation function Eq. (|57|) is real one 
obtains 



4 2 )(r)p = 2/(r) 



a-pa+ - ~ {o-+cr_,p} 



(64) 



so that up to second-order the corresponding master 
equation is indeed of the form of Eq. (|56|) with the expo- 
nential kernel function h(t) = 2f(t). However, in fourth 
order further terms appear which are not present in (|56l) : 



..(4) 



where 



it [e~ XT 



(1 - At) - e 



-2Arl 



(65) 



(66) 



As shown in [(| this implies in particular that if one trun- 
cates the expansion to first order in 70 the complete posi- 
tivity (and even the positivity) of the resulting dynamical 
map is violated for strong couplings in the Nakajima- 
Zwanzig case. On the contrary the second-order time- 
convolutionless master equation always guarantees com- 
plete positivity, as can be seen considering the second- 
order approximation for Eq. (|60j) given by 



7 (2) W =7o(l 



,-At 



(67) 



VI. CONCLUSIONS 



valid for a generic spectral density. Since the dynamical 
map giving the reduced system dynamics of the two-state 
system is known exactly, there is of course in principle 
no reason to resort to any kind of master equation in or- 
der to determine the dynamical behavior of the system. 
However, the present results lead to several important 
implications which are relevant for more realistic physical 
systems and their microscopic or phenomenological mod- 
elling, where analytical results cannot be obtained. In- 
deed for this model both time-convolutionless generator 
and Nakajima-Zwanzig kernel can be exactly expressed 
in terms of functions for which perturbative expansions 
are given, together with the exact solution for a reservoir 
with an exponential correlation function, corresponding 
to a Lorentzian spectral density. This allows for a de- 
tailed comparison of the two approaches expressing the 
dynamics in terms of a time-local and integrodifferential 
master equation respectively. It turns out that contrary 
to what is often expected the Nakajima-Zwanzig master 
equation is not simply obtained by convolution of the 
Lindblad operator appearing in the non-Markovian case 
with a suitable kernel. It actually has a different operator 
structure, emerging when considering higher perturba- 
tive orders. Furthermore the exact analytical result ob- 
tained for a Lorentzian spectral density shows the differ- 
ent convergence behavior of the two approaches. While 
the Nakajima-Zwanzig kernel is an analytic function of 
the coupling strength, providing a well-defined master 
equation at any time, the time-convolutionless generator 
breaks down at finite time in the strong coupling regime, 
thus failing to reproduce the asymptotic behavior. 



We have constructed the exact Nakajima-Zwanzig 
memory kernel for a specific model describing the de- 
cay of a two-level system into a reservoir of field modes 
which is initially in the vacuum state. The construction 
of the memory kernel is based on the analytical solution 
of the Schrddinger equation within the Hilbert space sec- 
tor describing states with zero or one excitation, and is 
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Appendix A 



In this Appendix we shall consider how to obtain a perturbative expansion for the function k\ (r) which according 
to Eqs. (|36]) and (j40l) - (j4T|) determines the memory kernel in the Nakajima-Zwanzig master equation, as a function of 
the two-point correlation function f(t) of the reservoir. To this end one considers the solution of Eq. (j2T)]) which is of 
the form Eq. flU with G (0) = 1 and G^ (t) explicitly given by 



G {2n) (t) = {-Y 



dU 



dU 



2i-l 



(Al) 



so that 



(i) = \G(t)\ 2 = 1 + 2K^T (-)" / *i / *a---/ *2nn/(*2i-i -*a 



JO JO 



J2(-Y l dt l dt 2 --- dt2nT[f(t2i-l-t2i) 
n=l ^ 70 ^ i=l 



'0 JO Jo 

Considering terms up to fourth order in the expansion parameter one has 

ft r ti 



z(t) = 1- 2K / dti / dhf(te) 
Jo Jo 

+2M [ dh [ dt 2 [ dt 3 f dt 4 f(h-t 2 )f(t 3 -t 4 ) 
Jo Jo Jo JO 



f dt x / dt 2 /(t 2 ) 

JO 



and denoting real and imaginary parts of f(t) as in Eq. fj42[) also 



z(t) = 1-2 f dti [ efta/i (t 2 ) 
Jo Jo 

+2 / dt x / dt 2 / dt 3 f dt 4 [fi (ti - t 2 ) A (*s - t 4 ) - /a (ti - h) h (*3 - *4)] 
Jo Jo Jo Jo 

/ dh [ dt 2 f 2 (t 2 

JO JO 



* ftl 

dh / dt 2 /i(t 2 ) 
Jo 



Introducing the functions 



hi(t) 



f dh f dt 2 fi(t 2 ) i = l,2 
Jo Jo 



one obtains for the Laplace transform of z{t) 

u — 2uf\(u) 2 r<£ a. i 



z(u) = 

u 



+ ^3 (-M U ) - /*(«)) + h\{u) + hl( U ) + 



k(u) = 2h{u) + - (h(u) + f 2 {u)) - u 2 (h\{u) + h 2 2 (u)) + 



and thanks to Eq. P5j) 



Using now the fact that the functions hi are equal to zero together with their derivatives at t = so that 

d^~ 



one has 



T ftl 



fc x (T) = 2/ x (r) + 2 / dtr/ dt 2 fi{h-t 2 )h(t 2 ) + 2 dh dt 2 f 2 (h-t 2 )f 2 {t 2 ) 



o Jo 

2 



T /.tl 



(*l) 



dhh(h) 



2/i(t) / dti / 1 dt 2 /i (ti - t 2 ) - 2/ 2 (r) f dh f 1 dt 2 f 2 (h ~ t 2 ) + 

JO JO JO JO 



o Jo 



T ftl 



We now exploit the identity 



r />i2 

dt 2 / dt 3 f{t 2 -t 3 )f{t 3 )- 
o Jo 



dt 2 fi (t 2 ) 



+ I dt 2 j 2 dt 3 f(T-t 3 )f(t 2 ) = 



o Jo 



9 



which can be checked noting that the function of t defined by the left-hand side of Eq. (|A10[) has vanishing derivative 
and is equal to zero for t — 0. We are thus left with 

k x {T) = 2/i (r) - 2 f dt 2 j 2 dt 3 [/1 (r - t 3 ) h (t 2 ) + /i(r)/i (t 2 - t 3 )} (All) 
Jo Jo 

-2 f dt 2 [' dt 3 \h (t - t 3 ) f 2 (t 2 ) + / 2 (r)/ 2 (f 2 - t 3 )] + • • • 
Jo Jo 

and thanks to the fact that real and imaginary parts of f(t) are even and odd respectively 

k x (t) = 2/j (t) - 2 / dt 2 / 2 di 3 [/1 (r - 1 3 ) /1 (-ta) - /a (r - i 3 ) /a (-fa) (A12) 
Jo Jo 

+/i(t)/i (t 3 - h) - h(j)h (*s - ta)] + •• • 
= 2/i (r) - 28 / dt 2 / 2 di 3 [/ (r - t 3 ) / (-i a ) + / (r) / (i 3 - fa)] + • • ■ 



Jo Jo 

Upon the change of variables i 2 — » i 2 — fi, f 3 — >• £ 3 — fi one has for the second and fourth order contribution to k 1 (r): 

fcf°(t-ti) = 2/i(t-ti) (A13) 

£^ 4) (t - ti) = -23? f dt 2 I* dt 3 [f (t - t 3 ) f (h -h) + f{t-t 1 )f (i 3 - t 2 )] . (A14) 
Jti Jti 

Appendix B 

Here we derive the contributions up to fourth order to the memory kernel Eq. (|36l) employing the standard Nakajima- 
Zwanzig projection operator technique. Since the initial state of system and bath is of the factorized form Eq. (|14p 
we can employ the standard projection operator 

Tw = Tr E [w)®p E , (Bl) 

where w is a state of system plus environment and pe denotes the vacuum state of the reservoir. This projection 
operator is the same used to obtain Eq. (fM)) and for it the initial state Eq. (fl"4")l is indeed an eigenoperator. Introducing 
further the superoperators 

C(t)p(t) = -»[#,(*),/,(<)] (B2) 

with H](t) as in Eq. ©, and 



£?(t,ti) = Texp QT dsQ£( 



(«) J (B3) 

where T denotes time ordering and Q = 1 — V, the Nakajima-Zwanzig memory kernel appearing in Eq. (|35[) is given 
by 

Knz(<-*i)p(*i) = "&B(A*)^(*.*i)aC(*i)p(*i)®PB)- ( B4 ) 
Noting that for this model 'PZI (fi) . . . £ (f 2n +i) "P = one has 

^nz(*-*i)p(*i) = TM>C(*)£(*i)/>(*i)®PiO (B5) 

+ / di 2 f 2 di 3 [Tr s (t 2 ) £ (t 3 ) £ (tx) p (ti) ® p B ) 

Jti Jti 

- Tr B (£(«)£ (t 2 ) PC (t 3 ) £ (fr) p (<i) ® + . . . 



Using Eq. © and Eq. P^|) one readily obtains 

= ( 

-/ {t - h) o-p {h) a+ + f {h -t)p (h) a+a- 



Tr E (C(t)C {h)p (h) ® p E ) = {-if [f (t - h) a+a-p (h) - / (f 1 - t) a^p (h) a+ (B6) 
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so that the second order contribution is given by 



(B7) 



+2/i (r) 



which due to Eq . (|46|) confirms the result Eq. (l52j) . Setting 

h(t,h,t 3 ,h)p(h) = Tr B {C(t)C (ta)jC(t 3 )r(ti)p(*i)(»pB} 

and 

h{t,h,h,tx)p{h) = r Tt B {C{t)C{t 2 )VC{h)C{t l )p{tx)®p E } 
a lengthy but straightforward calculation leads to the results 

/a (*, * 2 , *3, *i) p(tx) = / (t - *2) / (t 3 - *i) ct+^-P (tj) + f(t 2 -t)f (h - t 3 ) P (h) <r+a- 
+23? [/ (t - i a ) / (tj - t 3 )\ a+a-p (*i) cx+<r_ 
-4/i (i - i 2 ) A (<i - ts) ff-P (*i) <r+ 

and 

h(t,t 2 ,t 3 ,ti) p(h) = h (t,t 2 ,t 3 ,tx) p (tx) 

-23? [/ (t - t 3 ) f (h -t 2 ) + f(t- h) f (i 3 - t 2 )} a-p (ti) a+ 
+23? [/ (* - t 3 ) / (tx - ta) + / (t - t x ) / (t 3 - t 2 )] a+a_p (t x ) a+a_. 

One thus have for the fourth order expression 

Knz ( i - *i) P = -2B / *2 / 2 dt 3 [/ (* - * 3 ) / (*i - t 2 ) + / (t - ti) / (*a - t 2 )] 



(B8) 
(B9) 
(BIO) 



(Bll) 



(B12) 



which according to Eq. (1461) confirms Eq. ([543 
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